Instabilities in the two-dimensional cubic nonlinear Schrodinger equation 
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The two-dimensional cubic nonlinear Schrodinger equation (NLS) can be used as a model of phe- 
nomena in physical systems ranging from waves on deep water to pulses in optical fibers. In this 
paper, we establish that every one-dimensional traveling wave solution of NLS with trivial phase 
is unstable with respect to some infinitesimal perturbation with two-dimensional structure. If the 
coefficients of the linear dispersion terms have the same sign then the only unstable perturbations 
have transverse wavelength longer than a well-defined cut-ofi^. If the coefficients of the linear disper- 
sion terms have opposite signs, then there is no such cut-off and as the wavelength decreases, the 
maximum growth rate approaches a well-defined limit. 
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I. INTRODUCTION 

The two-dimensional cubic nonlinear Schrodinger 
equation (NLS) is given by 



(1) 




where tp = ip{x,y,t) is a complex- valued function, and 
a, j3 and 7 are real constants. Among many other situ- 
ations, NLS arises as an approximate model of the evo- 
lution of a nearly monochromatic wave of small ampli- 
tude in pulse propagation along optical fibers 
aP > 0, in gravity waves on deep water 0] 
aP < and in Langmuir waves in a plasma 
a(3 > 0. As a description of a superfluid 
is known as the Gross-Pitaevskii equation |^ ^ 
a(3 > 0. Sulem and Sulem ^ examine NLS in detail. 

NLS admits a large class of one-dimensional traveling 
wave solutions of the form 

ip{x, y, t) = (t){ax + by- st)e»At+»asa;+»6sy+»r,^ ^2) 

where is a real-valued function, s = s/{2aa'^ + 2/36^) 
and a, b, s, A and rj are realparameters. By making use 
of the symmetries of NLS ^ , all solutions of this form 
can be considered by studying the simplified form 



(3) 



where is a real-valued function, z ~ x — 2aKt, and n 
and A are real parameters. 
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If q;7 > 0, then NLS admits the following two solutions 
of the form 



0(z) = J2- fc cn(z,fc), with X ^ a{2k^ - 1 - k^) , (4) 



0(z) = y'2- dn(z,A:), with A = a(2 - fc^ - k^). (5) 
If a7 < 0, then NLS admits the following solution 



(j){z) = Y -2- k sn(z, k), with A = -a{l + k^ + K^). (6) 

Here k G [0, 1] is a free parameter known as the ellip- 
tic modulus and cn(-, k), dn{-,k), and sn(-, A:) are Jacobi 
elliptic functions. Byrd and Friedman provide a com- 
plete review of elliptic functions. If fc < 1, then each 
function 4i{z) is periodic. As fc ^ 1, the period of each 
increases without bound, and 4i{z) limits to an appropri- 
ate hyperbolic function, which we call a "solitary wave." 
These solutions, plus the "Stokes' wave" (plane wave) 



ipix,t) = Ae*' 



-i(qk — 7|A| )t 



(7) 
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comprise the entire class of bounded traveling wave solu- 
tions of NLS with trivial phase ■ Davey and Stewart- 
son show that a Stokes' wave is unstable unless either 
a/37 = Oj ov a(3 > and 07 < 0. In the remainder of this 
paper, we concentrate on Q), © and and on their 
instabilities. 

Zakharov and Rubenchik (ill establish that and 
with k — 1 are unstable with respect to long-wave trans- 
verse perturbations. Pelinovsky |13| reviews the stabil- 
ity of solitary wave solutions of NLS with af3 < and 
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a7 > 0, and presents an analytical expression for the 
growth rate of the instability near a cut-ofF. Extensive 
reviews of the stability of solitary wave solutions are given 
m The periodic problem has not been studied 

in as much detail, though Martin, Yuen and SafFman 
examine numerically the stability of the solution given in 
(O for a range of parameters. 

We present four main results in this paper. First, ev- 
ery one-dimensional traveling wave with trivial phase is 
unstable with respect to some infinitesimal perturbation 
with two-dimensional structure. For all choices of the 
parameters, there are unstable perturbations with long 
transverse wavelength. This generalizes the result of [ll| . 

Second, if a(3 > 0, then the only unstable pertur- 
bations have transverse wavelength longer than a well- 
defined cut-off. 

Third, if a/3 < 0, then there is no such cut-off. There 
are unstable perturbations with arbitrarily short wave- 
lengths in both transverse and longitudinal directions. 
These short wavelength instabilities seem to have been 
overlooked in previous analyses. 

Fourth, for a/3 < 0, the unstable perturbations with 
short wavelength have transverse wavenumbers that are 
confined to narrower and narrower intervals as the trans- 
verse wavenumber grows without bound. In these unsta- 
ble intervals, as the transverse wavenumber grows with- 
out bound, the maximum growth rate approaches a well- 
defined limit. As fc ^ 1, this limiting growth rate tends 
to zero if > 0, and to a finite non-zero limit if < 0. 



where p is a real constant, f2 is a complex constant, 
U and V are complex-valued functions and c.c. denotes 
complex conjugate. This leads to 

{aK^ + X)U - 3-f4>'^U + Pp^U - adlU = DV, (11a) 



{an^ + X)V - -fcf^V + Pp^V - adlV = -W. (lib) 

These are the central equations in this paper. We assume 
that U and V are periodic with the same period as (I). 
More general boundary conditions are discussed in \Vt\ . 
Instability occurs if (|ll|l admits a periodic solution with 
Re{Vl) < 0. Without loss of generality, for the remainder 
of this paper we assume k = by redefining A. 

In Sections IIIII and Hvl we examine (|lll) using small-p 
and large-p asymptotic analyses respectively. In Section 
Ivl we present results from a numerical study in which 
(|ll|l was solved for a wide range of p values. 



III. SMALL-p LIMIT 

Generalizing the work in |TH , we assume that for fixed 
k and for fixed small p, (|ll|l admits solutions of the form 



U ~ uo{z) + pui{z) + p'^U2{z) + 



(12a) 



II. STABILITY ANALYSIS 



V ^ vn{z) + pvi{z) + p^V2[z) + 



(12b) 



We consider perturbed solutions, -i/'p = V'p {x, y, t), with 
the following structure 



■0p = i(l)iz) + eu{x,y,t)+iev{x,y,t) + 0{e'^))e 



iK,x-{-iXt 



, (8) 



where u{x,y,t) and v(x,y,t) are real-valued functions, e 
is a small real parameter, z = a; — 2aKt, and (f){z)e^'^'^'^'^^* 
is one of the solutions presented in the previous section. 
Substituting ^ into 1^, linearizing and separating into 
real and imaginary parts gives 



^ p^oji + p^uj2 + 



(12c) 



where the Wj are complex constants and the Uj and Vj are 
complex- valued periodic functions with the same period 
as 4>. 

This assumption leads to the "neck" mode 



C/„(z,p) = 0(p), 



(13a) 



-(a/t^ -|- X)u + 3j(p^u + (3uyy + aux 



vt, 



(9a) 



K(^,p) = + O(p), 



(13b) 



— (aK^ -I- X)v + "f<j)^v + (3vyy + avxx = —ut- (9b) 

Without loss of generality, assume that u{x,y,t) and 
v{x, y, t) have the forms 



= -aPp'uJin + O(p'), 
and the "snake" mode 



(13c) 



u{x, y, t) = C/(z, p)e'py-^^ + c.c. 



(10a) 



C/,(z,p) = -f + 0(p), 
az 



(14a) 



i;(x,2/,i) =F(z,p)e''''«""* + c.c., 



(10b) 



y.(z,p) = 0(p), 



(14b) 
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FIG. 1: Plots of ujin and ujis versus k in (a) and (b) re- 
spectively. The solid line corresponds to (0J, the dotted line 
corresponds to l^l and the dashed line corresponds to 
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FIG. 2: Plots of IDmin /|a| versus k. The solid line cor- 
responds to ijlj, the dotted line corresponds to ||KJ and the 
dashed line corresponds to 



IV. LARGE-p LIMIT 

If a/3 > and p is chosen to be large enough to satisfy 

P'>5|^|, (15) 

then the two operators on the left side of Hill) have the 
same sign, so fi^ < 0. Therefore, there is no large-p 
instability if ajd > 0. 

If aP < and p is large, then one can show that there 
is no instability unless = 0(1). Therefore, we assume 

U r^Cl{pz) + p-^C2{^^z) + ■■■ , (16a) 



V ^U^Jiz) + p-^U^^z) + ■■■ , (16b) 



~ wi -I- p H , (16c) 



ap'^ = -l3p^ + v + 0{p-^), (16d) 

where j/ is a real constant, the Wj are complex constants, 
and the Q and are complex- valued periodic functions 
with the same period as (f). 

Substituting 1161) into Hll|l . one finds at leading order 

C/-Ciisin(/iz-f Zo) + O(p^^), (17a) 



= -a(3p^u;is + 0(p3), (14c) 

where LUin and wis are functions of the elliptic modulus 
of the unperturbed solution. Complicated but exact ex- 
pressions for ujin and ujis are derived in jT3|. The final 
results are presented in Fig. ^ where we plot ujin and 
Wis, the growth rates, versus A: for (|SJ| and 0. 

These plots establish that ujin < for Q and (|3Jl and 
Wis < for ®. Therefore, if a/? > 0, (Hj) and 10) are 
unstable with respect to long-wave transverse perturba- 
tions corresponding to the neck mode and © is unstable 
with respect to long-wave transverse perturbations cor- 
responding to the the snake mode. 

These plots also establish that ojis > for |0J and Q 
and ujin > for Therefore, if a/? < 0, Q and (0) 
are unstable with respect to the snake mode and © is 
unstable with respect to the neck mode. 

It follows from these results that a trivial-phase solu- 
tion of NLS is unstable to a growing neck mode if > 0, 
and to a growing snake mode if (3"f < 0. 



V ^^iism{pz + zo) + Oip-^), (17b) 

where Ciii Cii ^^'^ are constants. Requiring U and 
V to have the same period as (j){z) forces p to take on 
discrete values: p = 2TrN/L, where L is the period of 
(j){z), and TV is an integer. To satisfy (|16dp . iV > 1. At 
the next order in p, solutions are periodic only if 

wi = ±v/7/ - A + i^V^ - ^if - V, (18) 

where / is the Fourier coefficient (in sn\{pz + zq)) of 
4>^{z) sin(/iz -|- zq) and v is 0(1) but otherwise arbitrary. 
Minimizing the negative root in H18|) with respect to v 
leads to 

w„ii„ = -|7/|, (19) 

when V ~ {\— 2jf). Then (|16dp defines p„, the value of 
p at which the Nth unstable mode achieves its maximum 
growth rate: 

-PpI = a{27rN/Lf + 27/ - A. (20) 

We also find how for p can deviate from p„ before 
becomes negative 

Sp,^\^\f/(2\p\p,)^0{l/N), for 7V»1. (21) 
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Analytic expressions for the / corresponding to the 
solutions given in © and © are not known. But, 
in the large-/? limit, the Riemann-Lebesgue Lemma 
can be used to determine approximate expressions. As 
p oo, for the solution given in Q), 

^ ^ 2if(fc) I ^1 (^(''™(^^^^))) + " ^^^(^0 ' (^^^ 
for the solution given in ©, 

and for the solution given in JBJ, 

^ ^ W(fc) I ^ I (^^(^^ " ■ (24) 

In each of these expressions, am(-) gives the Jacobi ampli- 
tude and K{-) and E{-) are the complete elliptic integrals 
of the first and second kind respectively 0. 

Plots of — Wmm/|a|, a growth rate, versus k are given 
in Fig. 121 This argument establishes that all finite-period 
one-dimensional trivial phase solutions are unstable with 
respect to arbitrarily short-wavelength transverse pertur- 
bations if a/3 < 0, and that the growth rate of the insta- 
bility remains bounded as p — > cx). 

Note that as fc ^ 1, (j>{z) in (0 approaches a hyper- 
bolic tangent and the growth rate appr oaches that of the 
Stokes' wave with an amplitude of y^—2a/'j. This es- 
tablishes that there are an infinite number of of unstable 
branches if a/3 < and aj < 0. 

Also note that as k — + 1, 0(z) in both and Q 
approaches a hyperbolic secant, and the corresponding 
growth rate limits to zero. This establishes that there is 
no large-p instability in the solitary wave limit if > 0. 



The growth rates obtained from numerical simulations 
corresponding to © with k — -v/cTS and — a = /3 = 7 = l 
are included in Fig. O as dots. The line is obtained 
from the small-p results. The dashed curve is obtained 
from the large-p results with = 5. Each dotted 




FIG. 3: Plots of —Q versus p corresponding to with 
k — V 0.8 and ~a = (3 — ^ = 1. See text for a description. 
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FIG. 4: Plots of U and V versus z corresponding to (p = 
3.5, J7 = -0.99, = 5). 



V. MONODROMY 

The system of equations in Ijlll) is Hamiltonian in z, 
with periodic boundary conditions. The coordinates on 
the phase space are pi — dU/dz, p2 = —dV/dz, qi = U 
and q2 = V. The Hamiltonian is 

H = \{pI-pI) - + ~ il^'^i^M 

+ ^{X + PP^ -j4,^)ql + -q,q2. (25) 
2a a 

Such a Hamiltonian system necessarily has a mon- 
odromy structure with invariants 19]. We used this 
structure to identify the periodic solutions of by nu- 
merically integrating over one period of 0. 



curve corresponds to a different unstable mode. A plot 
of the spatial structure of the mode corresponding to 
(p = 3.5, n = -0.99, A^ = 5) is given in Fig. H 

Figure O demonstrates strong agreement between the 
numerical results and the small-p analysis when p is near 
zero. It also demonstrates agreement between the nu- 
merical results and the large-p analysis. 

Figure 01 demonstrates that the U and V obtained nu- 
merically are similar in form to the U and V obtained in 
the large-p analysis. 
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